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We study the suppression of electron localization due to the screening of disorder in a Hubbard- 
Anderson model. We focus on the change of the electron localization length at the Fermi level within 
a static picture, where interactions are absorbed into the redefinition of the random on-site energies. 
Two different approximations are presented, either one yielding a nonmonotonic dependence of the 
localization length on the interaction strength, with a pronounced maximum at an intermediate 
interaction strength. In spite of its simplicity, our approach is in good agreement with recent 
numerical results. 
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I. INTRODUCTION 



Understanding the interplay between disorder and 
electron-electron interactions remains one of the major 
challenges in modern condensed matter physics, experi- 
mental as well as theoretical. The research in this field 
has been stimulated by the possible metallic behavior 
in two-dimensional disordered interacting systemsii A 
metallic phase at zero temperature in two or less dimen- 
sions would be in contrast to the prediction of the scaling 
theory of Anderson localizationi^ The possible existence 
of a metallic phase, induced by interactions, is a long 
standing^ and still controversial problem, discussed by 
many authors during the last three decades 
One of the ideas proposed and discussed by several au- 
thors is that interactions lead to a partial screening of the 
random potential and, thus, reduce the effect of local- 
ization. In particular, the 2d disordered Hubbard model 
(the Hubbard- Anderson model) has been studied, mostly 
numerically, and it was demonstrated that repulsive in- 
teractions can have a delocalizing eSect^^^i^^^i^^^^^^ 
In this paper, we present an analytical study of the 
screening effect, focussing on the case of strong disor- 
der at zero temperature, when the Hubbard- Anderson 
model is in the regime of an Anderson insulator. Our 
approach is based on an exact treatment of the atomic 
limit, followed by "switching on" the intersite hopping t, 
under the assumption that the atomic-limit occupation 
numbers do not change. Let us emphasize that in our 
static approach, the interactions only change the original 
on-site energies. Therefore, we are left with a single- 
particle Anderson Hamiltonian, and the question is how 
the localized single-particle states may change due to the 
new, renormalized probability distribution of the on-site 
energies. In this sense, the approach is close, although 
not identical, to the Hartree-Fock treatment. A compar- 
ison of both methods will be presented. Although this 
approach is formulated for an arbitrary filling factor, it 
becomes inadequate close to half filling where magnetic 
effects dominate^ (such effects are not considered in our 



work). 

Furthermore, we will show that, for fixed disorder, the 
localization length ^ is a nonmonotonic function of the 
Hubbard interaction energy U , with a maximum for some 
intermediate value of U . This is because for strong 
interactions, the Mott-Hubbard physics of interaction- 
suppressed hopping dominates, leading to the forma- 
tion of two disorder broadened Hubbard bands with 
a reduced average density of states at the Fermi level 
and, as a consequence, increasing effectively the disor- 
der strength. In spite of the simplicity of the approach, 
the results are in good agreement with recent numerical 
studieSf^iii^i^iiSiii^ii^ in the appropriate range of parame- 
ters. 

The evaluation of the localization length within our 
present study is limited to one dimension. However, since 
the competition of screening on the one hand and Mott- 
Hubbard physics on the other hand operates in any di- 
mension, the nonmonotonic dependence of ^ on [/ should 
also hold in two and three dimensions, as argued below. 



II. ATOMIC-LIMIT APPROXIMATION 

In this paper, we consider the Hubbard- Anderson 
model, with on-site repulsion and on-site disorder, at zero 
temperature. The corresponding Hamiltonian is 



H = Ho + Hi 



kin 



H,. 



<i,j>.,(7 



(1) 



As usual, c|^(Cj^) denote fermion creation (destruction) 
operators of an electron at site i with spin a, n^^ = cl^c^^, 
t is the nearest-neighbor hopping amplitude, U is the 
on-site repulsion, ^ is the chemical potential, and 
are the on-site energies. The latter are assumed to be 
independent and uniformly distributed over the interval 
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FIG. 1: Site occupation in the atomic ground state with dou- 
bly occupied (black) , singly occupied (gray) and empty states 
(white), (a) weak interaction, (b) strong interaction for less 
than half filling (i.e., p < 1), and (c) strong interaction for 
more than half filling (p > 1). 



[— Y, y]i with the disorder parameter A. To focus on 
the screening effect in the case of strong locahzation, i.e., 
A ^ t, the interaction term will be absorbed into the on- 
site energies, yielding a renormalized distribution of the 
Ei. This results in an effective single-particle problem 
with a probability function pA{£i) which is derived as 
follows. 

In the atomic limit (< = 0), the ground state of the system 
can be solved exactly for an arbitrary filling factor p = 
^ , where Ne , N arc the numbers of electrons and lattice 
sites, respectively!^ The chemical potential jj, and the site- 
dependent occupation numbers {ni)o can be expressed as 
functions of p, A, and U: All sites with on-site energies 
below p—U are doubly occupied, all sites within [p — U, p] 
are singly occupied, and all other sites are empty (see Fig. 
[T]). Thus, the total occupation number of site i and the 
chemical potential are, respectively, 



{ni)Q 



andi^ 




p-U < El < p 

£i> P 



p <l, U < Ap 
p<l,U >Ap 
p>l,U<2A-Ap 
p>l,U>2A-Ap. 



(2) 



(3) 



The renormalized site energies depend on the site occu- 
pation numbers and can be read off from the poles of the 
(timc-ordcrcd) single-particle propagator. 



{ni,a)o{ni-a)o 



- p + U)-iO+ 

t)q) , (1 - {ni,a)o){'ni,^a)Q 



Lu — {si — p) — iO+ uj — {si 
UJ — {si — p) + iO+ 



p + U) + ^0^ 



(4) 



The first term in Eq. ^ corresponds to the doubly occu- 
pied sites showing that these on-site energies are shifted 
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FIG. 2: The renormalized on-site energy probability functions 
(solid line) for weak repulsion U : (a) atomic-limit approxima- 
tion, and (b) Hartree-Fock approximation. For comparison, 
also the original function is shown (dashed lines). 



by U . The next two terms correspond to singly occu- 
pied sites. In the absence of spin polarization, {ni_a-)Q — 
('^j,-o-}oj half of these on-site energies are again shifted by 
U whereas the other half remain unchanged. Finally, the 
last term corresponds to the unoccupied sites, whose en- 
ergies also remain unchanged. Combining Eqs. ([3]) and 
the rule for replacing the bare site energy Si by a 
renormalized one is 




each with 
prob. of i 



(5) 



For a weak to intermediate repulsion U these shifts lead 
to a raise of the lowest lying on-site energies towards 
the Fermi level p resulting in a renormalized probability 
function PAi^) with a reduced width and modified shape 
(see Fig. |2(a)| . 

In the atomic-limit approximation the Hamiltonian 
Eq. Il]) is replaced by the effective single-particle An- 
derson Hamiltonian 



H = 



t 



E 

<i,j>(J 



C „C, 



(6) 



with on-site energy probability function pa{£), where the 
two-particle interaction U enters only as a (screening) 
parameter in pA(e)- Note that, due to the asymmetry 
of pa{£), the average value of the renormalized on-site 
energies is 



J pA{£)£de 



(7) 
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In deriving Eq. ([6]), it was assumed that in the case of 
strong disorder, the occupation numbers in the atomic 
ground state, Eq. ([2]), are close to the occupation num- 
bers {rii) in the true ground state, with finite hopping 
amphtude t. This assumption is based on the results 
from the single-particle theory of localization,— where 
it is known that for strong disorder, an electron, once 
located at any site i, will stay at that site with high 
probability. More precisely, the change of the occupa- 
tion number of site z, 5{ni) = {rii) — {riij^, is of order 
t^/A. The same estimate holds if one considers hopping 
of a single electron on the background of the other elec- 
trons, which are assumed to be immobile. Furthermore, 
for strong disorder, a perturbative expansion in t, around 
the atomic ground state of Eq. ([T]), is possible.— Thus, 
with the same reasoning as in the noninteracting case, 
the leading self-energy corrections arc again only of sec- 
ond order. 

One measure of localization is the localization length ^ 
which governs the exponential decay of the single-particle 
wave functions ip{r) at distances far away from its local- 
ization center. It can be calculated from the probability 
for a transition from site i to site j as^^ 



lim 



log I \Gf A E)\ 



2 a;, 



(8) 



where G'^(i?) is the retarded propagator for a particle 
with energy E from site j to site i. In general, it is 
not possible to deduce ^ from the probability distribu- 
tion analytically, but in case of a one-dimensional lattice, 
there exists a relatively simple rclationi^ between ^ and 
{D{uj)), the disorder averaged density of states. 



CHE) = / {D{e)) log\E~e\ de. 



(9) 




(a) 




FIG. 3: (Color online) Localization length ^ at the Fermi level 
as a function of repulsion U and lattice filling p for A = 15: 
(a) atomic-limit approximation and (b) Hartree-Fock approx- 
imation. 



(Here and in the following, we choose units where t — 1 
and measure ^ in units of the lattice spacing.) 
In the strong disorder limit, {D{e)) can be replaced by 
Pa{£ + so that the inverse localization length at the 
Fermi level {E = 0) is given by 



PA{e) log\e - /i| de 



(10) 



A plot of f as a function of U and p, for fixed disorder 
strength A = 15, is shown in Fig. 3(a) It can be seen 
that for each given filling factor, the localization length 
exhibits a pronounced maximum. This maximum can be 
calculated to appear at 



A 

3" 
A 



1 -t- 3p(2 -p)-l 



p{2^p) + 0{p\2-pf) 



(11) 



The reason for this nonmonotonic behavior is simple: 
A weak to intermediate repulsion U changes the on-site 
energy distribution from a rectangular distribution to a 
narrower one by shifting low on-site energies toward the 
Fermi level (screening), see Fig. 



2(a) In contrast, 



very strong on-site repulsion enhances the localization 
by a large broadening of the probability density. There- 
fore, in between, there will be some value for which 
the screening is optimal and the localization length ac- 
quires a maximum. Such behavior is expected, since a 
strong repulsion effectively suppresses hopping processes 
and leads to an accumulation of spectral weight in the 
upper Hubbard band. 



In Fig. 4(a) the localization length ^ is shown as a func- 
tion of U for p ~ \ (quarter filling). A similar, non- 
monotonic behavior was also found in recent quantum 
Monte Carlo simulation o^'^'^'^i^^ and a most recent sta- 
tistical dynamical mean field theory evaluationi^ of the 
problem, where the conductivity and the inverse partic- 
ipation ratioi of a finite system were calculated, respec- 
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tively. Identifying a maximum of conductivity with a 
maximum of localization length, we find in all cases a 
good qualitative agreement with our results. Further- 
more, we even find a reasonable quantitative agreement 
with our results for the points of maximal delocalization, 
U^, Eq. (HI]), and Uf, defined in Eq. ^ below. Thus, 
there is a strong correlation between the degree of screen- 
ing and the conductivity: optimal screening corresponds 
to maximal conductivity. In addition, our analytical re- 
sults indicate that much of this physics of the nonmono- 
tonic behavior can already be understood on the level of 
a static screening approximation. 

Our result appears to be at odds with the statement^ 
that screening alone cannot account for the nonmono- 
tonic behavior of the conductivity. We will discuss this 
point in the next section. For that discussion, we need 
the variance of the renormalized on-site energies which is 
given by 
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{e)l 



(12) 



A2 
12 



pU + pU^ 
1 - 2p'U 
1-pU + pU^ 
l-2pU 



3U^ 
pU^ 

pU^ 



P<1 
P<1 
P> 1 
P> 1 



U <p 

U>P 
U <2-p 
U>2-p, 



with 

U 



U/A, p 



-- 3p(2 - p) , 

Hk-ip- 



p 

m 



3p(l - p) , 



Here, a\ has a minimum at 



A 



9/5-/5 



i+|i-p| 



3 P ^ 3 

otherwise. 



(13) 



(14) 



which can be seen in Fig. |4(b)| for p - 



i and A 



15. 



So far, our calculations were restricted to one di- 
mension because a generalization of relation ^ to two 
or three dimensions is, in general, not possible. However, 
in the limit of strong disorder, (E) can be calculated 
in good approximation by taking into account only the 
direct pathi^ from i to j. Therefore, Eq. ([5]) can be 
considered as a one-dimensional problem, leading again 
to Eqs. ([9]) and (fTO|) . respectively. Hence, we conjecture 
that the general result, i.e., the nonmonotonic behavior 
of localization, holds also in two and three dimensions. 



III. SITE-DEPENDENT HARTREE-FOCK 
APPROXIMATION 

The second model which will be discussed in this paper 
is the site-dependent Hartrce-Fock approximation]^ In 
this single-particle approximation, each site has a single 
renormalized energy level, given by 



U 



K)o, 



(15) 
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FIG. 4: (a) Localization length ^ at the Fermi level, normal- 
ized by its noninteracting value, as a function of the on-site 
repulsion (7 for A = 15 and p = \. (b) Variance of the 
renormalized on-site energy distribution, also normalized by 
its noninteracting value. The solid curves show the results 
for the atomic-limit approximation and the dashed curves are 
the Hartree-Fock results. 



which in Eq. ([T]) corresponds to the replacement, 



U 



Here, the absence of any kind of magnetization was as- 
sumed. Note that for a local, energy independent inter- 
action J7, the Hartree and the Hartree-Fock approxima- 
tions arc identical. The average on-site occupation {n^ ) 
was taken to be the one of the atomic ground state, (n^)o, 
Eq. (21), according to the assumption of a stable atomic 
configuration. 

As in the atomic-limit approximation, the shift of the 
occupied on-site energies leads to a renormalized distri- 
bution, with probability function pH(e) (Fig. 2(b)). The 



screening effect now is even more pronounced because 
all singly occupied states arc shifted by y yielding a 
smaller width and a stronger increase of the probability 
to find a state around the Fermi level. The resulting effec- 
tive single-particle Hamiltonian is again given by Eq. 
however, with the probability function p}i{e). 
The corresponding plots for the localization length in 
Hartree-Fock approximation are shown in Figs. |3(b)| and 
4(a) respectively. Again, a pronounced maximum of ^ 
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arises. In this case, it is found at the value 



A 



P(2-P), 



(17) 



which up to order O (p^ • (2 — p)^) coincides with the 
resuh from the atomic-Umit approximation Eq. (jlip . The 
increase of the locahzation length is considerably more 
pronounced in the Hartree-Fock approximation due to 
its narrower probability distribution and especially its 
larger probability density around the Fermi level. In both 
approaches, the effect of screening becomes stronger with 
decreasing disorder. However, for small values of A, the 
stability of the atomic ground state becomes doubtful. 
The average and the variance of the renormalized on-site 
energy distribution are, respectively, 



(e)H = 




with the same abbreviations as in Eq. 
is found at the value 



A/2 
A 



else. 



131). Its minimum 



(20) 




FIG. 5: Inverse localization length as a function of A for 
p = i and (a,) U — 7 in the atomic-limit approximation and 
(b) [/ = 14 in the Hartree-Fock approximation, respectively. 



As mentioned above, it was argued in Ref. that the pic- 
ture of screening would be too primitive to explain the 
nonmonotonic behavior and the evidence for a metallic 
state found in the conductivity simulations by varying 
the repulsion strength U. The argumentation was based 
on the observation, which the variance was a featureless, 
monotonically decreasing function of U around the tran- 
sition point. Our results show that this reasoning is gen- 
erally not conclusive. Although the static, single-particle 
treatment does not allow for the occurrence of a metallic 
state, we find a strong enhancement and nonmonotonic 
behavior of the localization length ^ as function of U, 
whereas the variance is also only a monotonically decreas- 
ing function around the point of maximal delocalization. 
For strong disorder and a distribution which is not char- 
acterized by a single parameter (like in the present case, 
cf. Fig. [5]), there is no simple relation between ^ and 
the variance a^. Especially for lower fillings, our results, 
Eqs. mi), dm) and respectively, show that 

the values of U for which the maximum of ^ and the 
minimum of cr^ occur, can be separated systematically. 
Moreover, we find that the atomic-limit approximation 
and the Hartree-Fock approximation do yield close val- 
ues of ^, although the variances can differ strongly (see 

Fig. a. 

In Refs. |8| and Il3l . the inverse participation ratio was 
calculated as a function of the disorder strength A, and 
a nonmonotonic behavior with evidence for a metallic 
state was found. It was argued there that the screening 



picture would necessarily predict a monotonic increase 
of the inverse participation ratio with increasing A, ex- 
cluding screening as a possible explanation. Our model, 
contrary to that statement, exhibits such nonmonotonic 
behavior as well, as shown in Fig. [51 where is plotted 
as a function of the disorder strength A, for some fixed 
values of p and U. This nonmonotonicity is caused by 
the crossover from the regime of disorder screening by 
interaction, A ^ t/, to the regime of interaction-reduced 
hopping, A <C C/, which is controlled by the parameter 
[//A. However, the exact position of the nonmonotonic- 
ity depends also on t/A, p, and the probability function. 



IV. CONCLUSION 

We examined the effect of static disorder screening 
by on-site repulsion in the one-dimensional Hubbard- 
Anderson model for strong disorder. We presented two 
different approximation schemes by absorbing the inter- 
actions into a redefinition of the single-particle on-site 
energies. In both approaches a renormalized probability 
distribution with an enhanced probability of finding site 
energies close to the Fermi level was obtained. We cal- 
culated the localization length at the Fermi energy for 
these single-particle problems and found a pronounced 
maximum of the localization length for some interme- 
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diate value of the repulsion strength. This can be un- 
derstood as a consequence of the fact that the increase 
of the localization length ^ for small U {U < U^) and 
the decrease of ^ for large U {U > U^) have different 
physical origins, namely disorder screening and reduced 
hopping, respectively. Similarly, a change of the "bare" 
disorder A, for fixed repulsion, resulted in a nontrivial, 
nonnionotic dependence of the localization length on the 
disorder strength. In contrast to the case of weak dis- 
order, we found no significant correlation between the 
variance of the effective on-site energy distribution and 
the localization length. Our results, especially in the case 
of the Hartree-Fock approximation, are in qualitative and 
to some degree even quantitative agreement with recent 
numerical studies. By our analytic approach, it was pos- 
sible to investigate the static screening effects separately 
from dynamical (inelastic) processes. 
We gave an argument that the same behavior should 
also be found for strong disorder in two and three 
dimensions 1^ In three dimensions there might be an 
interesting possibility of an interaction induced metal- 
insulator transition. Such a possibility is based on the 
assumption that our results, obtained in the strongly lo- 
calized regime, can be extrapolated up to the mobility 
edge. Under this assumption, a noninteracting Anderson- 



localized system, whose Fermi energy is sufficiently close 
to the mobility edge, will become metallic upon switching 
on interactions by shifting the mobility edge across the 
Fermi level. Furthermore, when the interactions exceed a 
certain strength, the system would reenter the insulating 
phase in analogy to the nonmonotonic behavior of the lo- 
calization length in one dimension.— (Let us emphasize 
that throughout this paper we do not discuss the case 
of half filling, with its characteristic Mott's physicsj ^^i^^ ) 
A similar, and experimentally more relevant, effect could 
also happen under a change of the bare disorder A, as 
suggested by Fig. [5l since the relevant dimensionless pa- 
rameter is U/A. 
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